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Abstract: I review the recent developments in the analysis of cosmic flow data, in particular, 
latest results of bulk flow measurements, comparison between redshift and peculiar velocity 
catalogs with emphasis on the measured value of the (= f2^ 6 /fe) parameter, and matter 
power spectrum estimates from galaxy peculiar velocity catalogs. Based on these develop- 
ments, one can argue that most of the previous discrepancies in the interpretation of cosmic 
flow data, maybe with the exception of bulk flow measurements on scales }t 100h~ 1 Mpc, have 
either been resolved or fairly understood. 



0.1 Introduction 

Within the gravitational instability (GI) framework for the growth of cosmic 
structures, the peculiar velocity field of galaxies and clusters provides a direct 
and reliable probe of the matter distribution, under the natural assumption that 
these objects are unbiased tracers of the large-scale, gravitationally induced, 
velocity field. The GI paradigm requires that the linear peculiar velocity, v(r) 
- defined as the deviation from Hubble expansion - and linear mass density 
contrast, S m (r), be related to one another according to the local (differential) 
relation, 



where VL m is the matter overdensity parameter. Note that the peculiar velocity 
field is determined by the distribution of the matter with all its components 
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or its global (integral) counterpart 
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especially the dominant dark matter component. 

In order to measure peculiar velocities of galaxies and clusters, observers use 
a variety of distance indicators. Generally, these indicators relate two quantities, 
one among which is distance dependent, e.g., galaxy luminosity, and the other is 
distance independent, e.g., galaxy rotational velocity. The best known examples 
of such indicators are the Tully-Fisher |24| and Faber-Jackson [jl3| relations; over 
the last decade or so these and many other types of distance indicators have 
been used to measure cosmological distances. The availability of a large number 
of galaxy peculiar velocity catalogs, some of them with few thousands objects, 
have turned cosmic flows to one of the main probes used to study the large scale 
structure in the nearby universe. Here I'll concentrate on the following three 
statistical measures of the velocity field: 

1. The Bulk Flow: This measure, defined as the average streaming motion 
within certain volume, is probably the easiest statistic to estimate from 
the observed radial component of peculiar velocities. At the Cosmic Mi- 
crowave Background radiation (CMB) restframe, the bulk motions are ex- 
pected to converge to zero with increasing volume. The rate of convergence 
depends on the fluctuations in the matter distribution on various scales, 
i.e., the large scale matter fluctuations power spectrum. This dependence 
on cosmological models has motivated several attempts to measure the 
dipolc component of the local peculiar velocity field and to determine 
the volume within which the streaming motion vanishes. As of yet, bulk 
flow measurements have produced conclusive and consistent results only 
on scales ^ 60/i _1 Mpc, but failed to do so on scales ^ 100/i _1 Mpc (for 
recent works see |, | | g || 13 H ) ■ 

2. The mass power spectrum: Equation^ suggests that one can estimate the 
bias free, O^ 2 weighted, matter power spectrum directly from the mea- 
sured peculiar velocities. To date, likelihood analysis based estimations of 
the matter power spectrum exist for the Mark III J35|, SFI |Q and EN- 
EAR catalogs. All of these measurements has consistently produced 
power spectra with amplitudes larger than those measured by other data 
sets, e.g., galaxy redshift surveys. 

3. The f3 parameter: Peculiar velocities enable a reconstruction of the large 
scale matter distribution independent of redshift surveys. Therefore, one 
can use Eqs. g and | to compare the matter and 3D velocity distributions 
deduced from the measured radial peculiar velocities to those obtained 
from redshift surveys. This comparison requires biasing model which spec- 
ifies how galaxies follow the underlying total matter distribution. On the 
scales of interest, it is usually assumed that 

5 g = bS m , (3) 

where 5 g is the galaxy observed density contrast and b is the linear bias 
parameter. The comparison is used to: 1) test the validity of Eqs. |l| and 
namely, the basic GI paradigm; 2) test the linear biasing model; 3) directly 
measure the value of Sl m (|| @, [u| ^ ^| and Q). Until recently, 
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comparisons using Eq. [I] have systematically yielded j3 values larger than 
those obtained from using Eq. |^. 

As mentioned above, many have attempted to estimate these statistical mea- 
sures of peculiar velocities over the years. Until recently, the results they ob- 
tained have often been inconsistent with each other and with estimations from 
other data, depending on the specific peculiar velocity catalog at hand and the 
analysis methods. It is generally accepted that the main reason for the in- 
consistencies lies in the problematic nature of the distance indicators used to 
determine the peculiar velocities. First, the distance measurements carry large 
random errors, including intrinsic scatter in the distance indicator and measure- 
ment errors, which grow in proportion to the distance from the observer and 
thus become severe at large distances. Further nontrivial errors are introduced 
by the nonuniform sampling of the galaxies that serve as velocity tracers. In 
particular, the Galactic disk obscures an appreciable fraction of the sky, creating 
a significant "zone of avoidance" of at least 40% of the sky. When translated to 
an underlying smoothed field, these errors give rise to severe systematic biases. 
In light of these difficulties, the inconclusive results have led many to question 
the reliability of the peculiar velocity datasets as cosmological probes. 

In this article, I review the recent developments in these three areas and show 
that significant improvements have occurred in the field in the past few years. 
In addition, an argument is put forward that those developements lead to alle- 
viating most of the inconsistencies and to understanding, at least qualitatively, 
the cause of the remaining outstanding ones. 



0.2 Bulk Flow 

The bulk flow motion direction and amplitude on different scales are the simplest 
quantities to measure from peculiar velocity data. They provide constraints on 
the power-spectrum of mass fluctuations. Theoretically, the mean square bulk 
velocity within a sphere of radius R, is given by, 

(v 2 (R)) = iff / P(k)W 2 (kR)dk, (4) 
Z7r Jo 

where P(k) is the mass fluctuation power spectrum and W(kR) is the Fourier 
transform of a top-hat window of radius R. In the upper panel of Figure [l| 
the rms expected bulk velocity, Vb rms (R) — (v 2 (R)), is plotted against R for 
a standard ACDM power spectrum with the dashed lines representing the la 
cosmic scatter deviations expected within the given volume. Note that the large 
cosmic scatter weakens the bulk flow statistic as a probe of cosmological models. 

On the observational side one can divide the currently available measure- 
ments into two main domains. The first is measurements of V;, within radius 
^ 60h~ 1 Mpc. Bulk flow values within this radius as measured from the SFI 
Mark III Q and the recently completed Shellflow [| and ENEAR || catalogs, 
lead to a roughly consistent picture even though some discrepancies still remain. 
For instance, the Mark III catalog yields a systematically larger amplitude of 
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Measured V b Direction 




Figure 1: Bulk Flow measurements. Upper panel: the symbols show the amplitude of 
the measured bulk flow (with its error) from the following surveys: Surface Brightness 
Fluctuations (SBF), SFI , ENEAR (EN), Shellflow (SF), Supernovae type la (SNIa), 
SMAC, EFAR, LP10, SCII and LP (see table for explanation) as a function of radius. 
The CMB dipole COBE measurement and bulk flow from the PSCz redshift catalog 
are also shown. The solid line shows the expected rms bulk velocity of a sphere of 
radius R for standard ACDM model; the dashed lines represent 1-cr cosmic scatter 
about the rms. Lower panel: the symbols show the direction of some of the measured 
bulk flow vectors, note that the catalogs that correspond to R ~ 60/i~ 1 Mpc have 
consistent directions while measurements that correspond to large distances do not. 
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Table 1: Recent Bulk Flow Measurements 



Survey 


v b 


R 


Comments 




(kms- 1 ) 


(/i^Mpc ) 




Tonry et al. |23[ 


350 


30 


SBF* 


Dekel et al. [fjjl 


370 


60 


M3* (TF 1 * + D n - a) 


Giovanelli et al. 


[p| 200 


60 


SFI (TF) 


Courteau et al. 


|l 70 


60 


Shellflow (TF) 


da Costa et al. \ 


f 220 


60 


ENEAR (Ai - cr) 


Riess et al. [20 


« 


100 


SNIa* 


Colless et al. (I 




- 100 


EFAR (FP°) 


Hudson et al. IT 


| 600 


140 


SMAC (FP) 


Lauer & Postman pj 700 


150 


LP (BCG*) 


Willick §§] 


700 


150 


LP10K (TF) 


Dale et al. |J 





- 150 


SCI/SCII (TF) 



* Surface Brightness Fluctuations method. * Mark III dataset. t Tully-Fisher Measure- 
ment. * Supernovae type la. Fundamental Plane measurement. * Brightest Cluster Galaxy 
method. 



the bulk motion ~ 370 kms -1 on scales ~ 60/i _1 Mpc as compared to values 
of < 220 kms -1 obtained from the other catalogs; this specific discrepancy is 
probably caused by the calibration problem the Mark III catalog is known to 
suffer from §, [TJ. 

Recently, results from the long awaited Surface Brightness Fluctuations 
(SBF) method have been first published p^ |. They give a bulk flow value of 
350km s -1 at 30/i -1 Mpc distance. The bulk flow measurement within this vol- 
ume, both in terms of amplitude and direction, is consistent with what is known 
about the structures just outside the sampled volume, e.g., the Great- Attractor 
the Perseus-Pisces superclusters. 

The second domain is measurements within R ^ 100/i -1 Mpc. The results 
here are far from being conclusive and various data sets lead to different bulk 
flow amplitudes and directions. While there is supporting evidence pointing 
towards convergence of the bulk flow on these large scales |j[ ||(], ||, other 
works [[l9], ^6[ argue for the existence of large amplitude ( £ 600 kms -1 ) 
streaming motions out to a depth as large as 150ft. _1 Mpc, ruling out that the 
Hubble flow has converged to the CMB frame at smaller distances. Given the 
far reaching implications that these large-scale motions would have on currently 
popular cosmological models it is clear that this issue is of great interest. It is 
important to point out, however, that the direction and the amplitude of the 
bulk motion detected on large scales by different authors do not agree with each 
other. Additional independent bulk flow estimates, e.g., from the PSCz redshift 
catalog and more recently from the dipole of the NVSS radio sources 0], are 
consistent with bulk flow convergence in the CMB restframe (see upper panel 
of Fig. 0). 

A summary of the bulk flow measurements is given in Fig. pi and Table [D. 
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Figure 2: Contour map of In- likelihood in the h — fl plane for ACDM models 
with 250 kms -1 thermal error component (added to account for the nonlinear 
environment in which early type galaxies tend to reside). The contours denote 
the most likely values within 1, 2 and 3a confidence levels. 



0.3 Power Spectrum Analysis 



The matter power spectrum has been measured from the Mark 111 | TM 3f| , SFI 
@ and ENEAR ]3lJ catalogs. With the exception of one study all the 
power spectrum estimations apply likelihood analysis which assumes that both 
the underlying velocity field and the errors are drawn from independent ran- 
dom Gaussian fields. The observed peculiar velocities constitute a multi- variant 
Gaussian data set, albeit the sparse and inhomogeneous sampling; this probabil- 
ity distribution function is reinterpreted as a likelihood function of the measured 
radial velocities given a model power spectrum. Maximizing the likelihood with 
respect to the model free parameters yields a best fit power spectrum. 

Together with its advantages (no explicit window function, weighting or 
smoothing the data and automatic underweighting of noisy, unreliable data) 
the likelihood method requires satisfying the following conditions: 1) explicit 
knowledge of the distribution function; 2) peculiar velocities are related to the 
over-densities through linear theory; 3) errors in the inferred distances constitute 
a Gaussian random field with scatter that scales linearly with distance. In 
addition, the likelihood analysis requires assuming a parametric functional form 
for the power spectrum. 

A typical result of the power spectrum parameters as estimated from, for 
example, the ENEAR data is presented in Fig. ||. This figure shows the like- 
lihood contour map in the TO — h plane, for the ACDM family of models 
with Harrison-Zeldovich spectrum (n — 1). In addition to the normal, distance 
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dependent, uncertainties the error matrix is assumed here to have a diagonal 
random contribution of 250 kms -1 that accounts for the nonlinear environ- 
ment within which the ENEAR early-type galaxies prefertionaly reside. The 
most probable parameters in this case (in the range Cl m < 1) are f2 m = 1 and 
h = 0.5. The elongated contours clearly indicate that neither tt m nor h are in- 
dependently well constrained. It is rather a degenerate combination of the two 
parameters, approximately Q m h x with x ~ 1 that is being determined tightly 
by the elongated ridge of high likelihood. 

As pointed out earlier, the power spectrum amplitude deduced from galaxy 
peculiar velocity data is considerably higher than those obtained from other 
types of data sets, usually favor the standard ACDM model (fi m = 0.3, A = 0.7 
and h = 0.65). This naturally raises the question of whether is the discrepancy 
driven by a yet undetected systematic effect in the data, or by something deeper? 

It is important to point out that the power spectrum estimation method is 
very sensitive to the assumed small scale power model (induced by noise and/or 
nonlinear evolution) which can add or suppress power. Recently, Hoffman and 
Zaroubi [^6| have carried out a detailed inspection of the goodness-of-fit of the 
Mark III, SFI, and ENEAR best-fit power spectra by using Principal Component 
Analysis (PCA) approach and found that in neither of the three cases do the 
assumed theoretical power spectrum and/or error model give an acceptable 
fit. In a subsequent study, Silberman et al. [p2j have shown that this misfit is 
probably caused by ignoring small scale power in the analysis; thereby casting 
a large shadow of doubt about the applicability of the error models or the strict 
validity of linear dynamic, assumed so far. 



0.4 The Value of /? 

Galaxy peculiar velocities and their redshift space positions have been used to 
estimate the value of (3 = Q® n e /b, under the hypotheses of linear theory and 
linear biasing. These analyses have been typically carried out using two alter- 
native strategies. In the so-called density-density comparisons a 3-D velocity 
field and a self-consistent mass density field are derived from the observed ra- 
dial velocities and compared to the galaxy density field measured from large 
redshift surveys, under the assumption of linear bias (a combination of Eqs. [l] 
& ||). The typical example here is the comparison of the density field recon- 
structed applying the POTENT method g, |l) to the MARK III catalog of 
galaxy peculiar velocities |^7| with the IRAS 1.2 Jy redshift catalog density 
field pl| . The various applications of density-density comparisons to a number 
of datasets have persistently led to large estimates of (3, consistent with unity 
pl| . The alternative approach is that of the velocity- velocity analysis. In this 
second case the observed galaxy distribution is used to infer a mass density field 
from which peculiar velocities arc obtained and compared to the observed ones 
(a combination of Eqs. || & ||). The velocity- velocity methods have been applied 
to most catalogs presently available yielding systematically lower values of /3, in 
the range of 0.4 — 0.6. 

Both density-density and velocity- velocity methods have been carefully tested 
using mock catalogs extracted from N-body simulations. They were shown to 
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Figure 3: /3 from density-density comparison. Left Panel: The 1, 2 and 3 a likelihood 
contours in the 13 — AS plane. Right Panel: A scatter plot showing the SEcat UMV 
reconstructed G12 density vs. the PSCz G12 density. In this comparison, one out of 
ten from all the data points that have reconstruction error < 0.2 was randomly picked. 
The size of the symbols is inversely proportional to their errors. 



provide an unbiased estimate of the /? parameter. Yet, when applied to the same 
datasets the discrepancy in the [3 estimates turned out to be significantly larger 
than the expected errors. Velocity-velocity comparisons are generally regarded 
as more robust as they require manipulation of redshift catalog data as opposed 
to the density-density comparison which manipulates the noisier and sparser 
peculiar velocity data. In any case both approaches are quite complicated and 
it is hard to understand how systematic errors can arise and propagate through 
the analysis. It is, therefore, likely that these systematic effects do influence the 
P parameter estimation. 

Since all the previous results have been obtained with methods designed to 
conclusively carry out either velocity- velocity comparison (ITF M] & VELMOD 
|H ||)) or density-density comparison (POTENT it is very difficult to 

know whether the source of the discrepancy is data or methodology driven. Re- 
cently, a new linear approach, called unbiased minimal variance (UMV) method 
has been proposed |gfj. With the UMV method, both velocity- velocity and 
density-density comparisons can be carried out within the same methodological 
framework. 

Similar to the Wiener Filter the UMV estimator is derived by re- 

quiring the linear minimal variance solution given the data and an assumed 
prior model specifying the underlying field covariance matrix. However, unlike 
the Wiener filter, the minimization is carried out with the added constraint of 
an unbiased reconstructed mean field. In the context of reconstruction from 
peculiar velocity data, the UMV algorithm could be regarded as a compromise 
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Figure 4: /3 from velocity-velocity comparison. Left Panel: The 1,2, and 3 a likelihood 
contours in the f3 — AH a plane from a comparison of the SEcat G12-smoothed radial 
velocities at the locations of the SEcat data points vs. the PSCz G12-smoothed radial 
velocities at the same locations comparison. Right Panel: A scatter plot of the G12- 
smoothed radial velocities, in this comparison all data points have been included. The 
size of the symbols is inversely proportional to the errors 



between the POTENT algorithm, which assumes no regularization but might 
be unstable to the inversion problem of deconvolving highly noisy data, and 
the Wiener filter algorithm, which takes into account the correlation between 
the data points and therefore stabilizes the inversion, but constitutes a biased 
estimator of the underlying field (for detailed discussion see |$(], |3^] ) . 

The UMV approach has been applied to the SFI, ENEAR, Mark III and 
to the newly compiled SEcat catalog (which is a combination of the SFI and 
ENEAR catalogs) |3^ ]. The reconstructed fields are compared with those pre- 
dicted from the IRAS PSCz galaxy redshift survey to constrain the value of (3. 
For example, the analysis of the SEcat catalog for the first time leads to consis- 
tent (3 values from the comparison of the density and the velocity fields yielding 
= 0.57lo.i3 an d P = 0-51 ± 0.06, respectively. The (3 values obtained from 
comparing the PSCz data to each of the four galaxy peculiar velocity catalogs 
are consistent. 

Figure || shows the estimated value of (3 from the PSCz and the SEcat 
density-density comparison. The fields where smoothed with 12/i _1 Mpc Gaus- 
sian kernel; only points with small errors, estimated from Monte Carlo mock 
SEcat realizations, where used in the comparison. The left panel shows a con- 
tour plot in the plane of f3 and A5 which corresponds to zero-point density 
offset. The 1, 2 and 3-cr level certainty contours are shown. These contours 
are drawn assuming that the error in the density at each individual point is 
independent, this assumption is obviously wrong. Therefore, in the (3 results 
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Table 2: /? Measurements 



Method 


Compared data 


& 


S — S Comparison 


POTENT @ 
UMV|| 


Mark III us. IRAS 1.2Jy 
SEcat vs. PSCz 


0.89 ±0.12 

u -°'-0.13 


v — v Comparison 


VELM 
VELM 
ITF |1 
ITF |7| 
UMV 


od ||,| 

OD §] 

D 

a 


| Mark III vs. IRAS 1.2Jy 
SFI vs. PSCz 
Mark III vs. IRAS 1.2 Jy 
SFI us. IRAS 1.2Jy 
SEcat vs. PSCz 


0.50 ±007 
0.42 ±0.07 

?t 

0.6 ±0.1 
0.51 ±0.06 



t Inconsistent flow fields (probably due to problematic calibration of Mark III) 

quoted in the previous paragraph the uncertainty has been estimated by adding, 
in quadrature, to the likelihood analysis 1-a uncertainties the error in recover- 
ing (3 as estimated from the Monte Carlo simulations. The right panel shows a 
scatter plot of densities with their best fit linear model. Figure || shows results 
similar to the previous figure, but for the reconstructed radial velocity field. 
The comparison here is somewhat different in the sense that it is carried out 
between the reconstructed 12/i _1 Mpc smoothed radial velocities at the location 
of the SEcat data points. For more details on the comparison see J?2| , 

Summary of the (3 parameter values obtained from various velocity-velocity 
and density-density comparisons is shown in table |^. 

Although, the results obtained by the UMV analysis strengthen the case 
for lower (3 values, by no account do they render the results of the Mark III- 
POTENT analysis invalid. However, they call for re-examination of the later 
bearing in mind the following issues. The Mark II catalog, as shown for example 
by Davis et al. [[l0| and more recently by Courteau et al. B, suffers from sys- 
tematic calibration error that would cause a systematic error in the estimation of 
(3. However, this error is not expected to affect the density-density comparison 
as to overestimate the value of (3 by more than a factor of two. An application 
of the UMV method to the Mark III catalog |3^] shows that the obtained values 
of (3 are much lower than those obtained by applying the POTENT analysis pl| 
but still they are somewhat higher than those obtained from the SEcat cata- 
logs by 0.1 — 0.2. Moreover, the velocity-velcotiy analysis (like VELMOD) yield 
consistent values of (3 when applied to Mark III and SFI datasets (|, ||^, |29|. 
Based on these arguments, one could speculate that the most likely explanation 
to the previous inconsistency is a result of a collusion between the systematic 
errors in Mark III and the noise sensitivity of the POTENT reconstruction. 

0.5 Summary &; Discussion 

To evaluate the outstanding issues in the peculiar velocity measurements in gen- 
eral terms one should differentiate between two types of problems. The first is 
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the inconsistencies among results obtained from various peculiar velocity data 
sets and occasionally even within those obtained from the same data set. The 
second, is the disagreement with measurements from other types of data sets. 
The former is obviously more severe as it implies unresolved systematics in the 
peculiar velocity datasets themselves. Measuring by this yard stick, the con- 
flicting results obtained from the bulk flow and the (3 parameter measurements 
are more serious than the higher amplitude found with the power spectrum 
measurements . 

It is reassuring that within the local universe (R £ 60ft.~ 1 Mpc) all the recent 
bulk flow measurements, especially from the SFI, Shellflow and ENEAR catalogs 
do agree very well with each other both in terms of amplitude and direction. 
Those three marginally agree with the Mark III data, which gives a somewhat 
higher value of Vf,. Given that the Shellflow sample have clearly shown that 
the Mark III data set is somewhat miscalibrated, this disagreement is hardly an 
issue. 

Unfortunately however, the picture on larger scales (R }t 100/i _1 Mpc) is 
very different and the disagreement among the various measurements is yet to 
be resolved. While the bulk flow obtained from the SNIa §(}, SCI/SCII § 
and EFAR J0| samples clearly points toward convergence of the CMB dipole; 
the SMACgTj, LP (T| and LP10 |§ surveys find a bulk flow amplitude of 
~ 600— 700km s _1 . Although the former three measurements are consistent with 
the one obtained from the PSCz redshift catalog data, the NVSS radio sources 
and with theoretical prejudice, by no means the later three have been refuted. 
It is worth pointing out however, that the LP bulk flow, while comparable in 
amplitude, disagrees with the SMAC and LP10 flow. Furthermore, almost all of 
the deep peculiar velocity surveys have small number of objects and therefore 
probably prone to the pitfalls of small number statistics. 

The peculiar velocity measurements have systematically led to mass power 
spectrum amplitudes higher than those obtained from other types of data. In 
light of the overwhelming evidence pointing towards lower amplitude power 
spectrum, a special effort should be made to rule out any inherent bias in the 
prior assumptions made in the peculiar velocity based power spectrum measure- 
ments. In fact, two recent studies ^2) strongly suggest a generic problem 
with the theoretical framework used to estimate the mass power spectrum from 
the Mark III, SFI, and ENEAR velocity surveys. Possible sources of this prob- 
lem lie with non- linear dynamical effects and/or oversimplified treatment of 
errors p2f . 

The irreconcilability of the [3 estimation from density-density and velocity- 
velocity comparisons has been one of the major outstanding issues in the cosmic 
flows field of study. A new technique, the UMV, have enabled carrying out, for 
the first time, both comparisons in the same framework. The results obtained 
by applying this technique to several galaxy peculiar velocity catalogs yield 
low values of the parameter (~ 0.5 — 0.6), a result consistent with those ob- 
tained from the previous velocity-velocity comparisons and from the analysis of 
redshift surveys. These latest results clearly strengthen the case for low (3 val- 
ues, in agreement with those obtained by the previous velocity-velocity studies. 
A qualitative inspection into the reason of the discrepancy in the estimation 
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of (3 between the UMV and POTENT density-density comparisons allows one 
to speculate that the most likely explanation is a collusion between both the 
systematic errors in the Mark III data and the noise effects on the POTENT 
algorithm, that somehow conspired to produce these high /3-values. 

In light of the developments presented in this review, it is argued that most 
of the outstanding issues in the large scale peculiar velocity field of study, maybe 
with the exception of the large scale velocity field, have been either resolved or 
understood (at least qualitatively) and we finally have a consistent cosmological 
model emerging from the study of cosmic flows. The experience gained during 
the convergence towards the current status, despite the crooked path it took, 
will be invaluable when dealing with the large future datasets. 
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